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ABSTRACT 

The quantum phenomena of spectral flow which has been observed in laboratory 
superfluids, such as '^He-B, controls the drift velocity of proton type II superconductor 
vortices in the liquid core of a neutron star and so determines the rate at which 
magnetic flux can be expelled from the core to the crust. In the earliest and most 
active phases of the anomalous X-ray pulsars and soft-gamma repeaters, the rates 
are low and consistent with a large fraction of the active crustal flux not linking the 
core. If normal neutrons are present in an appreciable core matter-density interval, 
the spectral flow force limits flux expulsion in cases of rapid spin-down, such as in the 
Crab pulsar or in the propeller phase of binary systems. 
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1 INTRODUCTION 

It is now well-established that evolution and decay of the 
internal magnetic field, rather than rotational energy, is the 
dominant source of the persistent and burst X-ray luminosi- 
ties in certain hot young isolated neutron stars. These are 
observed as the anomalous X-ray pulsars and soft-gamma 
repeaters, referred to as magnetars (for reviews see Harding 
& Lai 2006; Woods & Thompson 2006; Pons, Link, Miralles 
& Geppert 2007). In principle, a large part of the internal 
magnetic flux could be contained in the liquid core of the 
star at densities less than ~ 2po, where po = 2.8 x 10^* g 
cm~^ is the nuclear density. In this interval, many calcula- 
tions of pairing and of the equation of state (Baldo, Burgio 
& Schulze 1998; Heiselberg & Hjorth- Jensen 2000) show that 
the core consists of neutrons, protons, relativistic electrons 
and possibly muons, and that the neutrons and protons are 
superfluids, pairing in the F2 and ^So states, respec- 

tively. The proton energy gaps found in these calculations 
indicate a type II superconductor (Baym, Pethick & Pines 
1969) However, Baldo & Schulze (2007) have since shown 
that the inclusion of more reflnements, including three-body 
forces, reduces the energy gap to values consistent with type 
I superconductivity for certain proton densities. It is also 
the case that other authors have more recently considered 
the possibility of type I superconductivity with various mo- 
tivations (Sedrakian, Sedrakian & Zharkov 1997; Buckley, 
Metlitski & Zhitnitsky 2004; Sedrakian 2005; Alford, Good 
& Reddy 2005; Charbonneau & Zhitnitsky 2007; Alford & 
Good 2008). But it remains a reasonable assumption that 
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type II superconductivity exists for at least part of the den- 
sity interval we are considering here. 

In the type II system, the magnetic flux is quantized in 
units of cf>o = hc/2e — 2.07 x 10~^ G cm^ and is confined 
to the cores of proton vortices (see, for example, Tinkham 
1996). Steady-state motion of the proton vortices is deter- 
mined by balance of the Magnus force, electromagnetic inter- 
actions, external forces including buoyancy forces (Muslimov 
& Tsygan 1985; Jones 2006) and interaction with vortices of 
the rotating neutron superfluid (Muslimov & Tsygan 1985; 
Sauls 1989; Jones 1991; Ruderman, Zhu & Chen 1998), and 
forces derived from interactions between the quasiparticles 
of the system. The distinct populations of superfluid quasi- 
particles concerned are those localized in the cores of proton 
vortices and those delocalized in the continuum. Previous 
vortex velocity calculations (Jones 2006) assumed both the 
internal temperatures T < 10* K of typical radio pulsars 
and negligibly small neutron and proton quasiparticle num- 
ber densities. Under these conditions, proton vortices and 
therefore the internal flux distribution could be quickly ex- 
pelled to the core-crust boundary. 

The low-temperature assumption is not appropriate for 
magnetars whose X-ray luminosities in the active phase, 
~ 10'^ — 10** yr after formation, are thought to be a con- 
sequence of plastic flow of the crust (with its frozen-in fleld) 
under very high magnetic stresses (Harding & Lai 2006; 
Woods & Thompson 2006). This dissipative process must 
maintain a large input of thermal energy to the inner crust 
and core of the star so that internal temperatures are likely 
to be much higher than those inferred for the isolated radio 
pulsars. Neglect of continuum superfluid quasiparticles is not 
justified at these higher temperatures unless it can be es- 
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tablished that the neutron energy gap is largo at aU relevant 
wavenumbers. Although early calculations gave such values 
(Amundsen & 0stgaard 1985), much work (see the review 
by Heiselberg & Hjorth-Jensen 2000) finds a small angle- 
averaged '^P2 — F2 neutron energy gap, A„ ~ 0.1 McV, at 
the wavenumbers concerned. More recently, Khodel, Khodel 
& Clark (2001) also obtained similar small values for cer- 
tain neutron-neutron interactions. But there is also the pos- 
sibility of pion condensation at a matter density ~ 2po. Its 
effect has been investigated by Khodel ct al (2004) and is 
to increase the gap in the neutron single-particle spectrum. 
Nevertheless, we regard the existence of a small neutron su- 
pcrfluid gap, or possibly even a normal Fermi liquid, as be- 
ing not improbable for some wavermmbcr bands within the 
density interval considered in this paper. 

A new development in condensed-matter physics, the 
phenomenon of spectral fiow described by Volovik (1996; 
see also the review of Kopnin 2002), has greatly clarified and 
extended our understanding of the forces acting on moving 
vortices and has been confirmed by measurements on the 
^He-B phase (Bevan et al 1997). The present paper applies 
these ideas to show that forces generated by spectral flow 
limit the velocities of proton vortices in the outer core of hot 
young neutron stars to an extent that is consistent with a 
large fraction of the active flux not linking the core. 



2 SPECTRAL FLOW 

The fermion zero-mode excitations concerned are those 
quasiparticles localized within an isolated vortex core at 
energies below Ap (referred to the chemical potential as 
zero). Approximate eigenvalues (Caroli & Matricon 1965; do 
Gennes 1966) in the limit pp^ 3> 1, where ^ is the superfluid 
coherence length and pp the Fermi wavenumber, are £'p^, 
where /it is the (half-integer) angular momentum component 
parallel with the vortex axis z, and 



Epfj, = —p,TwJo, 



■Km*^^p± 2EfpPi_ 



Ti^Pf 



< Ap. 



(1) 



(This expression neglects interaction with the core magnetic 
flux; see Bardeen et al 1969). The proton bare and effective 
masses are nip and m* and its Fermi energy is Epp- The 
components of the proton wave- vector parallel with and per- 
pendicular to the vortex axis are p and Values of n are 
large at the temperatures of Table 1, so that localized quasi- 
particles interacting with an external heat-bath (consisting, 
in this case, of a distribution of continuum quasiparticles) 
can be treated by a Boltzmann equation (see Kopnin 2002) 
in which the quantum number /x is replaced by a continuous 
variable. In this semi-classical model, the vortex structure is 
that of a cylindrical hole in the superfluid order parameter 
and quasiparticle states within this hole arc normal particle 
and hole orbits in the form of chords limited at each end 
by Andreev reflection (Andreev 1964), the specific process 
of reflection that occurs at a normal-superfluid interface. 

The perturbation of the Boltzmann distribution func- 
tion caused by vortex motion consists of a displacement of 
the Fermi sphere in the plane perpendicular to z which also 
contains the vortex velocity vl. It relaxes through binary 
collisions with either localized (vortex-core) quasiparticles 



or with continuum quasiparticles. Thus the relaxation rate 
is given by the transport relaxation times, t^^ = t^I+Tc'^, 
but only the latter process transfers momentum between 
the vortex and the continuum. The existence of these two 
distinct relaxation processes moans that, within the semi- 
classical model, the localized quEisiparticles mediating the 
force between the continuum heat-bath with bulk velocity 
vt and a vortox of velocity v_l have a steady-state bulk 
velocity v;oc = 'vlt/tioc + •vtt/tc- This lags behind vl in 
general, and in the limits Tc or tioc 00 tends, as expected, 
to Vi or vr. Thus spectral flow (Volovik 1996) depends on 
the velocity difference 



V = Vl - Vioc 



(vl - Vt)t/Tc. 



(2) 



In Andreev reflection, a particle (hole) is retroreflected as a 
hole (particle), so maintaining the orientation of the chord. 
It has angular momentum n = [(r — vt) x p_l] • z, where 
r — vt is its coordinate in the vortex rest frame. The retrore- 
flcction is not exact owing to the motion of the superfluid 
so that the chord precesses about the vortex axis with an- 
gular frequency wo in a direction opposite to that of the 
circulating superfluid (Stone 1996). Motion of the vortex in- 
duces a finite rate of change of orbit angular momentum, 
given by the classical variable /x = — [v x px] • z, that drives 
spectral flow, which is the movement of quasiparticles both 
to and from the negative energy continuum Epp, < along 
the continuous dispersion relation given by equation (1). For 
— /u, ~ 10 or more, we rely on the correspondence principle in 
assuming that classical-orbit quasiparticles give an adequate 
description of spectral flow. In the case of the physical dis- 
crete spectrum, with half-integer /u, the process of spectral 
flow relics on the existence of sufficient level-broadening to 
allow quasiparticles to move along the dispersion relation by 
hopping from state to state. We refer to Stone (1996) and 
Kopnin (2002) for details of the relationship between the 
continuum and discrete-/! cases. Solution of the Boltzmann 
equation, followed by integration of the rate of change of 
quasiparticle momenta over fj, and px, gives the force act- 
ing on the vortex, including the Magnus and spectral flow 
components. We do not give further details here because 
the review by Kopnin (2002; pp. 1654 - 1660) contains a 
complete and very transparent account. 

On this basis, the force-balance equation per unit length 
of vortex is. 



fB + fv + Fm + Fi + Fsf = 



(3) 



in which fs,y are the external forces caused by buoyancy 
and interaction with neutron vortices. The lordanskii force 
(Thouless, Ao & Niu 1996; Sonin 1997) is unimportant here 
and has been neglected in equation (3). The Magnus force 
is expressed in terms of the proton electromagnetic current 
density Jp obtained from the entrainment coefficients (An- 
dreev & Bashkin 1975; Alpar, Langer & Sauls 1984; Jones 
1991), vl and the superfluid bulk velocities. 



Fm = -Jp X 00. 



(4) 



A similar expression gives F;, the electromagnetic interac- 
tion with electrons and muons, in terms of the lepton current 
density J;. The spectral flow force has components parallel 
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with and perpendicular to vl — vt, 
1 T 



-c 



1 + wgr^ Tc 



z X (vi - Vt), 



1^0 T 



Vt) 



(Volovik 1996; Kopnin 2002), with 



G ■ 



nhpp 



tanh 



(5) 



(6) 



(7) 



where = knT. The total proton density is Pp = Pp + pp, 
but at the temperatures considered here, the normal compo- 
nent Pp is negligibly small compared with the superfluid Pp. 
The simple expressions given by equations (5) and (6) are 
valid at temperatures such that loottc » 1, which we shall 
find to be always satisfied here. The definition of relaxation 
times follows Kopnin and Lopatin (1997). 

At this point, it is interesting to note that frictional 
forces on vortices in astrophysical superfluids have been es- 
timated previously from the scattering of continuum quasi- 
particles by core quasiparticlcs, with neglect of spectral flow. 
Following this procedure, Sedrakian (1998) found an expres- 
sion for the frictional force on neutron vortices arising from 
the scattering of continuum proton quasiparticles. In the 
present context, the order of magnitude of such a force would 
be 



(8) 



This is smaller than the spectral flow component, equation 
(6), by a factor of the order of ksT/Epp, where Epp is the 
proton Fermi energy, and so is negligible by comparison. It 
follows that the frictional force produced by the scattering 
of electrons or muons is also negligible because the transport 
rela^xation time for this process is several orders of magni- 
tude longer than the strong-interaction values of Tc- 



3 APPLICATION TO PROTON VORTICES 

Spectral flow is a remarkable quantum phenomenon whose 
recognition has enabled the understanding of laboratory 
measurements of the force acting on moving ^He-B vortices 
(see, for example, Bevan et al 1997). Its predictions are 
transparent in the semi-classical modelling of what would 
otherwise be an extremely cumbersome problem. The ratio 
of vortex radius to average interparticle spacing is given by 
Pf^ which, for ^He-B, has values in the interval 10^ — 10^. 
Typical values in the proton superfluid arc perhaps an or- 
der of magnitude smaller. For example, from the energy gap 
obtained by Baldo & Schulze (2007; the curve in Fig. 2 of 
that paper containing all corrections) we find pp^ = 7.5 at 
po = 0.16 fm~^, but this increases rapidly thereafter as a 
function of matter density. Use of the Boltzmann equation 
does require that the number of quasiparticles be large, and 
this condition is satisfied if we accept that the effective num- 
ber is that contained within a length of vortex at least an 
order of magnitude greater than ^. The circumstances exist- 
ing in the core of a neutron star also differ from those, for 
example, of Bevan et al (1997) in relation to the problem of 
defining the heat bath and its velocity. Apart from collective 
excitations, the heat bath contains the lepton-photon system 



and continuum quasiparticle excitations of the neutron and 
proton superfluids. It is also possible that the neutrons may 
exist as a normal Fermi liquid within limited intervals of 
density. All these components interact electromagnetically 
or strongly with the nuclei embedded in the solid crust of 
the star and, in the absence of vortex motion, would be sta- 
tionary with respect to it. But any component more strongly 
coupled with moving vortices than with the crust would be 
expected to flow with them. The neutrons are the most im- 
portant component and in their case we shall follow the ar- 
gument given by Kopnin (2002; Sect. 9 of that paper) in re- 
lation to the analogous problem in ^He-B and assume that 
the continuum quasiparticle rest frame is stationary with 
respect to the crust provided their mean free path for scat- 
tering by other continuum quasiparticles is shorter than the 
mean free path for scattering by localized proton quasipar- 
ticles. In a steady state of vortex motion, the force fs -|- fv 
is transferred to the neutrons and must be balanced by a 
neutron chemical potential gradient. But this is very small, 
of the order of 10"^ - 10"^ eV cm-\ 

Equation (4) has been expressed as shown to emphasize 
that Fm -|- Fi = because screening of the electromagnetic 
current in the whole volume of the neutron-star core must be 
presumed for a type II superconductor (Jones 1991, 2006) 
except for the microscopic circulating currents of the vor- 
tices. The distribution of leptonic current density is present 
throughout core and crust prior to the superconducting tran- 
sition and in non-superconducting regions, can change only 
on the same long time-scales as the magnetic flux density. 
Thus the superconducting transition must be to a state of 
finite supercurrent density such as to satisfy the condition 
Jp -|- Ji =0 everywhere in the core. In effect, the lepton- 
proton system is an insulator and the moving-vortex induc- 
tion field induces no large-scale current density. But Jp is 
confined to the core of the star, unlike J;, and the return cur- 
rents that complete the circuit flow as sheets on the surface 
of the spherical superconducting volume and exclude the 
magnetic flux present in adjacent non-superconducting re- 
gions. In having volume screening, neutron stars differ from 
laboratory superconductors for which screening is primarily 
a surface phenomenon (see Tinkham 1996). 

In the low temperature limit, ljot 3> 1, only a small 
longitudinal component of the spectral flow force, given by 
equation (6), remains. Owing to the cancellation of the Mag- 
nus force, equation (3) reduces to fB+fv+Fsfw = 0, showing 
that movement of proton vortices, and hence expulsion of the 
field, occurs easily under the influence of buoyancy forces or 
of interaction with outward-moving neutron vortices as the 
rotation of the star slows (Jones 2006). At higher temper- 
atures, the existence of spectral flow radically changes this 
behaviour. 

It is convenient to take account of both longitudinal and 
transverse components of the spectral flow force by obtain- 
ing an expression for the component of v_l — vt parallel with 
fs + fi'. Provided the relaxation times are such that equa- 
tions (5) and (6) are valid, this component is independent 
of T and is. 



(VL - Vt)|| 



OJoTc 

~~C~ 



(fs+fv). 



(9) 



Neutrons are the most important heat-bath component. The 
relaxation time Tc for interaction with superfluid neutrons. 
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wavevcctor k, is not well known owing to its exponential de- 
pendence on the •^P2 —^F2 energy gap and our obligatory ne- 
glect of the possible existence of nodes in this function and of 
its k-space anisotropy. Given these limitations, we make the 
most elementary possible approximation and assume that 
the neutron-proton G-matrix is isotropic and identical with 
the transition matrix for free-particle neutron-proton scat- 
tering at an equivalent neutron laboratory energy. We obtain 



n / 



(10) 



in which (7„p is chosen at a laboratory energy equal to the 
sum of the Fermi energies, Epn + Epp (for this, the amgle- 
averaged centre-of-mass energy in the medium equals that 
for the free-particle scattering). The function W \s W = 1 
in the case of a normal Fermi system, but for an isotropic 
superfluid is given by, 



1 

W 



T77 = dx 



dy 



1 -I- e-^ 1 + 



(11) 



in which = x'^ + (/3A„)'^, with an identical defini- 
tion of Y. The numerical values of W obtained in the in- 
terval appropriate for Table 1 axe satisfactorily fitted by 
W^ = 0.8(/3A„)-iexp(/3A„). 

We adopt effective masses m* — O.Smp and mj^ = 
0.9m„ (Zuo, Bombaci & Lombardo 1999) in evaluating equa- 
tions (1) and (10) and for reference purposes, assume a pro- 
ton fraction of 0.03 at po (Heiselberg & Hjorth-Jensen 2000) 
and thus a Fermi energy Epp = 7.4 MeV. The factor relat- 
ing the dissipative force in equation (9) to the size of the 
Magnus force is woTc. The Table 1 evaluations of this are 
for p± = pf and for Ap = 0.062 MeV, selected as approx- 
imately the minimum value consistent with type II super- 
conductivity at this proton density, which gives a reference 
value LOci — 1.23 x 10^* rad s~^. At the temperatures of the 
Table, the condition u)qttc S> 1 is satisfied because tioc has 
the same order of magnitude as the normal- neutron Tc- The 
question of the extent to which the neutron quasipaxticle rest 
frame moves with the vortices rather than being stationary 
with respect to the crust can be answered, for the temper- 
atures and An of the Table, by noting that anp and the 
proton-proton cross section app at the equivalent laboratory 
energy equal to 2Efp differ only by a factor of order unity. 
(See Yao et al 2006; pp. 339-340. For example, app « 150 
mb at 20 McV, whereas a^p ~ 100 mb at 70 MeV.) The 
Kopnin condition is approximately satisfied provided the ra- 
tio of the neutron continuum quasiparticle density to the 
spatially-averaged density of localized proton quasiparticles 
((X ApB^^ cxp(— /3A„)) exceeds unity. Dragging of the heat- 
bath rest frame is negligible for normal neutrons. In the case 
of superfluid neutrons at the lower temperatures and smaller 
A„ , it is significant only for spatially-averaged magnetic flux 
densities B ~ 10^* G, but is much reduced for Ap larger than 
the minimum. 

The effective component of the vortex drift velocity is 
given by equation (9) with C = 1.69 x 10^° erg s cm~^ at 
matter density po, neglecting the hyperbolic function. It is 
only weakly-dcpcndcnt on proton-system parameters other 
than Ap. Values for the case of stronger proton pairing can 
be found simply by scaling from the Table 1 reference value 
of ujo using the new values of Ap and Epp- In particular, it 



Table 1. Values of ojqTc are given in the second column, as a func- 
tion of temperature, for interaction of localized (vortex-core) pro- 
ton quasiparicles with normal neutrons. The remaining columns 
give this quantity in the case of superfluid neutrons with an as- 
sumed isotropic A„. Vacant spaces are for temperatures close to 
or above the critical temperature, or for temperatures so low that 
LUoTc is too large to be of interest. The assumed value of uio is the 
reference value, 1.23 x 10^** rad s~^. 



T 

(10«K) 



An 

(MeV) 





0.00 


0.04 


0.08 


0.12 


0.16 


0.5 


420 










1.0 


105 


1910 








1.5 


47 


271 


2971 






2.0 


26 


94 


477 


3240 




2.5 


17 




150 


641 


3060 


3.0 


12 




68 


212 


750 


3.5 


8.6 




37 


94 


265 


4.0 


6.6 




23 


50 


119 



is to be emphasized that there is no exponential dependence 

on the value of the proton energy gap. 

The buoyancy force /s is of the order of Hci<I>o/8'itR ~ 
1.0 dyne cm^^ (Muslimov & Tsygan 1985; Jones 2006), 
where Hd is the superconductor lower critical field and R is 
the neutron-star radius. This term can arise because the nor- 
mal to superconducting transition is so fast that the super- 
conductor is in a steady, but not necessarily complete equi- 
librium state. One aspect of this is that, even for B < Hci, 
the magnetic fiux cannot be expelled from the proton super- 
conductor (Baym, Pethick & Pines 1969). Even if fs were 
negligible, the term fv would exist in the case of superfluid 
neutrons, caused by interaction between neutron and proton 
vortices (Muslimov & Tsygan 1985; Sauls 1989; Jones 1991; 
Ruderman, Zhu & Chen 1998). It is limited by a critical 
value fy at which neutron vortices cut through the proton- 
vortex lattice, whose magnitude is at most of the same order 
as the buoyancy term (Jones 1991). If the outward radial 
velocity of neutron vortices during spin-down exceeds the 
velocity given by substituting fs -I- fy in equation (9), they 
will move by cutting through the proton vortex lattice in 
those cases where a relative sliding motion is not possible. 

General predictions are, of course, impossible to make 
owing to the exponential dependence of spectral flow on A„ . 
Together with Ap, this function of density also controls the 
core neutrino cmissivity and hence the internal temperature 
of isolated neutron stars less than ~ 10^ yr in age (see for 
example. Fig. 5 of Yakovlev & Pethick 2004) which lack any 
substantial source of internal dissipation. It is quite possible 
that values of A„ larger than those assumed in Table 1 ex- 
ist in certain density intervals, but it is the regions of small 
An that have the greater effect on average drift velocities. 
The choice of ^P2 — ^ F2 neutron energy gaps in the Table 
is consistent with recent calculations (Heiselberg fc Hjorth- 
Jensen 2000; Dean & Hjorth-Jensen 2003). At kp ~ 2.0 
fm"'^ and above, higher partial waves make little contribu- 
tion to attractive forces and it is even possible that there are 
signiflcant intervals in which the neutrons are not superfluid. 
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4 CONCLUSIONS 

It is of some interest that a new quantum phenomenon rec- 
ognized in relation to laboratory studies of vortices in ''He 
should have some bearing on neutron-stax magnetism. Sub- 
ject to the reservations described in Sect. 3, the conclusions 
of this paper that follow from the spectral flow phenomenon 
are as follows. The presence of charged or neutral baryons as 
a normal Fermi system within some interval of core-matter 
density leads to wqTc values of the order of those in the 
seond column of Table 1. These are some orders of magni- 
tude smaller than the loqTc ~ 5 x 10^ needed for velocities 
~ 3 X 10~^ cm s~^ that give significant magnetic flux trans- 
fer from core to crust within the 10"' yr magnetax phase. 
The absence of substantial movement of flux from core to 
crust in this time interval would be consistent with a field 
configuration in which the distribution of the active fiux be- 
lieved to be responsible for magnetar phenomena is confined 
to the crust at or close to the time of neutron star format 
tion, its ability to move and evolve being unconstrained by 
linkage with the core. This is also true, even in the absence 
of normal neutrons, for the A„ of Table 1 because the higher 
temperatures there may well be maintained by internal dis- 
sipation during the active magnetar phase. Comparison of 
the continuum and localized quasiparticle densities shows 
that there can be dragging of the heat-bath frame vt by 
vl, much reducing spectral fiow, particularly for very high 
fields, B ~ 10^® G, and for Ap near the minimum consistent 
with type II superconductivity, but our broad conclusions 
are unchanged. 

There are other circumstances in which spectral fiow 
may be of relevance, such as the case in which the density- 
dependence of neutron-pairing allows both normal and su- 
perfluid neutrons within appreciable intervals of matter den- 
sity. The values of looTc given in the second column of Table 
1 for modest internal temperatures T « 10* K are small 
enough to constrain the proton vortex drift velocities that 
are possible in cases of rapid spin-down such as the Crab 
pulsar, and the propeller phase of binary systems. The out- 
ward movement of neutron vortices in the supcrfluid region 
then occurs by sliding relative to the proton vortices or by 
intersecting them, and is unable to force the outward move- 
ment of magnetic flux. 
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